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ABSTRACT 


Travel-time  and  magnitude-yield  calibration  of  Northern  Eurasia,  which  is  largely  aseismic,  can  be  improved  by 
using  the  large  chemical  and  peaceful  nuclear  explosion  (PNE)  seismic  data  sets  acquired  by  the  Soviet  deep  seismic 
sounding  (DSS)  program.  Eleven  major  data  sets  of  this  program  have  recently  been  digitized  and  have  become 
available  to  nuclear  test  monitoring  research.  Although  the  PNEs  have  been  already  used  in  several  monitoring 
studies,  the  chemical-explosions  data  sets  are  still  incompletely  explored.  In  the  present  study,  we  use  inversion  and 
numerical  modeling  to  analyze  the  coda  amplitude  parameters  and  P-wave  travel  times  from  numerous  DSS 
explosions,  which  could  lead  to  a  detailed  empirical  mapping  of  seismic  properties  within  northern  Eurasia. 

For  proper  mapping  of  coda  parameters,  it  is  important  to  select  a  parameterization  that  reflects  the  fundamental 
physics  and  is  stable  in  respect  to  unconstrained  factors.  We  investigated  two  types  of  such  parameterizations, 
namely  the  frequency-dependent  coda  attenuation  Qcoda(f)  and  frequency-independent  attenuation  with  geometric 
spreading  (%QC0 da)-  By  utilizing  a  newly  developed,  parallelized  version  of  the  one-dimensional  (1-D)  reflectivity 
program,  we  were  able  to  model  within  the  0.2  to  10  Hz  frequency  band  and  to  800  sec  times,  which  nearly  covers 
the  complete  DSS  time-frequency  range.  Several  complexly  layered  velocity  models  related  to  the  study  area  were 
tested. 

The  numerical  modeling  showed  that  for  a  fixed  intrinsic  S-wave  Q  of  the  crust,  the  resulting  Qcoda(f)  exhibited  a 
strong  positive  frequency  dependence.  This  spurious  frequency  dependence  is  due  to  a  disregard  of  ray  bending  and 
other  geometrical  propagation  effects  and  not  related  to  model  rheology.  Notably,  similar  strong  Qcoda(f)  ©cy0  8_  1 0 
dependencies  were  observed  in  many  PNE  coda  records  and  also  in  other  areas.  By  contrast,  coda  amplitude 
inversion  using  the  geometric- spreading  model  results  in  a  stable  spreading  parameter,  y~  -0.008  and  a  frequency- 
independent  gCoda  varying  from  ~2Q  for  Q  =  300  to  ~1.2 Q  when  Q  =  1000.  The  variations  in  Qcoda(Q )  as  well  as  the 
fact  that  2coda>  Q  are  consistent  with  a  part  of  seismic  coda  energy  propagating  through  the  mantle.  Tests  for 
variations  of  the  crustal  velocity  structure  suggest  that  the  Qcoda(Q)  dependence  is  relatively  stable,  and  therefore  it 
can  be  used  to  invert  for  the  intrinsic/scattering  attenuation  of  the  crust  from  coda  measurements. 

We  have  assembled  a  database  of  coda  records  from  selected  PNEs  and  chemical  explosions  and  performed  their 
inversion  and  preliminary  interpretation  based  on  the  (%QC0 da)  model.  As  in  the  modeling,  parameter  ^is  relatively 
stable,  whereas  Qcoda  shows  variations  that  can  be  associated  with  tectonic  structures.  Crustal  Q  parameters  are 
further  mapped  from  the  inverted  modeled  Qcoda(Q )  dependence. 

In  order  to  further  expand  the  DSS  database,  we  have  obtained  earthquake  recordings  from  several  DSS  PNE 
profiles.  In  the  ongoing  work,  we  are  utilizing  these  data  for  coda  magnitude  analysis  and  derivation  of  empirical 
magnitude  (apparent  coda  source  spectra)  relationships.  Additionally,  we  present  an  update  of  our  first-arrival 
travel-time  calibration  model  for  Northern  Eurasia  using  the  PNE  as  well  as  selected  chemical-explosion  records. 
These  efforts  should  ultimately  lead  to  a  unified  three-dimensional/two-dimensional  (3-D/2-D)  calibration  model  of 
Northern  Eurasia,  including  travel  times,  stable  coda  magnitudes,  and  magnitude-yield  relations.  This  model  will 
bridge  the  gap  between  Europe,  central  Asia,  the  European  Arctic,  and  eastern  Siberia  and  provide  additional 
calibration  information  for  further  studies  in  these  areas. 
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OBJECTIVES 

Many  seismic  monitoring  studies  are  based  on  properties  of  the  seismic  coda  from  earthquakes  and  explosions  and, 
in  particular,  on  its  amplitude  decay  and  attenuation.  Short-period  coda  amplitude  decays  lead  to  the  most  stable 
estimates  of  source  magnitudes  (Phillips,  2004),  and  measured  coda  attenuation  (Q)  is  often  related  to  crustal 
properties  that  are  difficult  to  measure  otherwise  (Aki  and  Chouet,  1975).  High-frequency  (>5  Hz)  codas  from  PNEs 
in  the  east  European  platform  also  led  several  researchers  (Ryberg  et  al.,  1995;  Enderle  et  ah,  1997;  Morozov  et  ah, 
1998a,b;  Ryberg  and  Wenzel,  1999;  Morozov,  2001)  to  propose  sharply  contrasting  models  of  small-scale 
heterogeneity  of  the  uppermost  mantle. 

Despite  its  broad  use  in  seismological  practice,  the  causes  of  variations  of  coda  Q  with  frequency  are  still  not  well 
understood.  In  many  studies,  strong  frequency  dependencies  of  coda  Q ,  as  well  as  body-  and  surface-wave  Q ,  were 
documented.  This  dependence  is  commonly  presented  in  a  power-law  form: 

Q{f)=Q0[yfj ,  CD 


where  f0  is  some  reference  frequency  often  taken  to  equal  1  Hz  (Aki  and  Chouet,  1975).  Both  Q0  and  exponent  77  are 
considered  constant  within  the  frequency  band  of  interest.  Note  that  Q  typically  increases  with  frequency  (77  >  0), 
and  large  values  of  77  ~  0.5-1  are  not  uncommon  in  Lg  Q  and  Lg  coda  Q  studies  (e.g.,  Nuttli,  1973;  Hasegawa,  1985; 
Der  et  ah,  1986;  Campillo,  1987,  1990;  Benz  et  ah,  1997;  Frankel  et  ah,  1990;  Mitchell  and  Cong,  1998;  Mitchell  et 
ah,  1997,  1998;  McNamara  et  ah,  1996;  McNamara,  2000;  and  Erickson  et  ah,  2004).  Visco-elastic  rheological 
models  were  proposed  to  explain  observations  of  frequency-dependent  Q  in  Earth  materials  (Liu  et  ah,  1976). 
Correlation  to  tectonics  in  North  America  suggested  that  active  tectonic  regions  are  generally  characterized  by  low 
Qo  and  high  77,  while  stable  cratons  are  generally  characterized  by  higher  Q0  and  lower  77  (Erickson  et  ah,  2004). 
However,  in  our  PNE  studies  in  the  Siberian  Craton,  extremely  high  values  of  77  ~  1  were  found  in  combination  with 
very  high  Q0  (Morozov  et  ah,  2006).  For  77  ~  1,  the  coda  amplitude  decay  becomes  in  fact  frequency-independent, 
and  this  suggests  a  hypothesis  that  high  apparent  77  could  in  fact  indicate  low  attenuation  and  significant  geometric 
spreading.  An  argument  in  favor  of  such  an  interpretation  will  be  given  below. 

In  coda  measurements,  its  Q(f)  value  is  derived  from  the  following  model  of  time-frequency  dependent  coda 
amplitudes: 


A(t,f)=A0(f)ex  p 


Q(f) 


(2) 


where  t  is  the  time  after  the  primary  arrival  and  A0(f)  is  the  coda  spectrum  at  7=0.  Note  that  the  key  assumption  of 
model  (2)  is  that  the  geometrical  spreading  of  coda  waves  (°c  in  Aki  and  Chouet,  1975)  is  compensated  by  the 
scattering  volume  increasing  as  the  time  t  increases.  Similarly,  note  that  in  surface  wave  studies  (e.g.,  Erickson  et  ah, 
2004),  frequency-dependent  Q  can  be  measured  only  by  time-domain  techniques,  which  also  rely  on  compensation 
of  geometrical  spreading  (Nuttli,  1973). 

In  this  report,  we  show  that  observations  of  frequency-dependent  coda  Q  may  be  related  not  to  crustal  rheology  but 
to  the  chosen  parameterization  Eq.  (2).  Compensation  of  geometrical  spreading  in  Eq.  (2)  should  be  relatively 
complete  for  local  earthquake  codas  involving  straight  scattered  rays  and  uniform  scattering  volumes  (Aki  and 
Chouet,  1975).  However,  at  crustal-scale  and  regional  distances,  scattering  tends  to  concentrate  near  the  surface,  and 
scattered  waves  consist  of  complex  assemblages  of  diving,  refracted,  reflected,  guided,  and  surface  modes.  Based  on 
various  assumptions  about  these  modes  and  stochastic  properties  of  the  crust,  several  coda  models  can  be  used 
(e.g.,  Aki  and  Chouet,  1975;  Sato  and  Fehler,  1998).  In  practice,  however,  with  virtually  any  coda  model,  additional 
ad  hoc  distance  and  spectral  corrections  are  required  (Mayeda  and  Walter,  1996;  Phillips  et  ah,  2004).  Here,  we 
introduce  such  a  correction  by  noting  that  an  uncompensated  (most  likely,  undercompensated)  geometrical 
spreading  G(t )  should  remain  in  the  coda.  Because  it  is  only  practical  to  seek  a  two-parameter  coda  description 
[similar  to  ( Q0 ,rj)  in  Eq.  (1)],  we  abandon  the  frequency  dependence  of  2coda  in  relation  (2)  and  approximate  the 
geometrical  spreading  as  exp(-jf),  yielding: 
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4*>/)=4>(/)ex  p 


f  1 

— 

r+ 

t 

\  J~  coda  y 

(3) 


where  the  geometric  spreading  (relaxation)  parameter  /is  measured  in  units  of  1 1  time.  A  similar  expression  for 
frequency  dependence  of  surface  wave  attenuation  was  given  by  Dainty  (1981).  Note  that  the  extreme  case  of  7J  =  1 
[Qocfm  Eq.  (2)]  simply  corresponds  to  geometrical  spreading  with  low  attenuation  (2coda=°°)  in  relation  (3). 

Our  objective  below  is  to  show  that  the  (Q^rj)  description  of  frequency-dependent  coda  attenuation  (and  maybe  also 
in  other  types  of  bulk  crustal  Q)  can  be  effectively  and  usefully  replaced  with  (%  gcoda)-  The  principal  advantages  of 
such  a  replacement  are  three-fold:  1)  formula  (3)  does  not  rely  on  assumptions  about  geometrical  spreading  but 
empirically  parameterizes  it  instead;  2)  geometrical  spreading  and  Q  are  explicitly  separated;  and  3)  0coda  becomes 
frequency-independent,  in  agreement  with  simple  Newtonian  solid-state  mechanics,  at  least  within  the  presently 
available  data  resolution.  The  somewhat  puzzling  increase  of  Q  with  frequency  suggesting  depletion  of  scatterers  at 
shorter  scale  lengths  becomes  explained  by  geometrical  spreading  due  to  diving  (among  other)  scattered  waves.  In 
addition,  as  our  analysis  of  PNE  profiles  shows  (see  example  below  and  Part  II  of  this  paper),  parameter  ^appears  to 
be  relatively  stable  across  large  areas,  whereas  Qcoda  shows  strong  variability  between  different  tectonic  structures. 

Below,  we  analyze  the  ( Q0,r j)  and  (y  Qco da)  parameterizations  of  coda  attenuation  using  numerical  examples.  We 
generate  four  finite-difference  models  of  PNE  coda  wavefields  resulting  from  realistic  crustal  and  upper-mantle 
velocity  models  and  test  them  for  variations  of  crustal  velocities  and  Q.  As  expected,  the  ( Q0 ,  //)  parameterization 
leads  to  strong  spurious  Qcoda(f)  dependencies.  Inversion  of  this  Qcoda(f)  for  the  true  crustal  Q  could  represent  a  non¬ 
trivial  task.  By  contrast,  in  the  (y  Qco da)  model,  the  values  of  ^can  be  predicted  from  the  crustal  structure,  and  0coda 
can  be  uniquely  related  to  Q  of  the  crust.  Repeated  for  several  types  of  velocity  structures  and  crustal  velocities,  this 
modeling  allows  us  to  construct  “calibration”  dependences  Q(QC0 da)  from  numerical  models  and  use  them  to  invert 
the  observed  coda  data  for  crustal  properties  in  future  research.  In  addition,  we  derive  an  approximate  mapping  of  (y 
Gcoda)  parameters  into  ( Q0 ,  //)  and  discuss  the  potential  causes  of  observing  low  or  high  values  of  77.  It  appears  that 
rather  than  7j  or  Q0,  as  an  attribute  of  the  velocity  structure,  ^could  be  a  better  factor  differentiating  between  stable 
(cold)  and  active  (hot)  tectonic  zones. 

Overall,  it  appears  that  parameterization  (y  Qco da)  is  not  only  safer  theoretically,  but  also  provides  the  basis  for 
inverting  coda  amplitudes  for  bulk  crustal  S-wave  attenuation  and  for  tying  them  to  the  tectonics.  Because  of  its 
explicit  separation  of  geometrical  and  attenuation  effects,  the  approach  is  also  likely  to  help  in  building  more  robust 
and  transportable  coda  magnitudes  and  0-based  seismic  regionalization. 


RESEARCH  ACCOMPLISHED 

We  examine  the  effective  frequency  dependence  of  coda  attenuation  by  numerically  modeling  the  teleseismic  P 
wave  coda  resulting  from  scattering  of  seismic  energy  from  near-surface  crustal  heterogeneities.  The  use  of  near¬ 
surface  crustal  scatterers  simplifies  our  computations,  and  it  is  also  justified  by  the  increased  heterogeneity  of  the 
shallow  crust  and  supported  by  several  array  observations  (e.g.,  Greenfield,  1971;  Dainty,  1985  and  1990;  Bannister 
et  al.,  1990;  Gupta  et  al.,  1991)  and  previous  coda  modeling  (Dainty  and  Schultz,  1995).  The  choice  of  the  P-wave 
coda  is  due  to  several  reasons  (Morozov  and  Smithson,  2000):  1)  S-  and  Lg  phases  are  less  pronounced  in  the  PNE 
records;  2)  their  codas  lie  on  top  of  the  P-wave  coda,  complicating  the  observation  and  requiring  their 
decomposition;  and  3)  the  Lg  phase  is  significantly  more  difficult  to  model  accurately.  Similarly,  the  Rayleigh  wave 
( Rg )  decays  quickly  in  near-surface  sediments,  and  it  is  not  viewed  as  a  significant  factor  in  coda  formation  (Dainty, 
1985).  Rg  is  not  seen  beyond  -200  km  in  PNE  records  and  is  also  difficult  to  reproduce  numerically;  consequently  it 
is  not  examined  here.  The  use  of  P-wave  coda  of  arrivals  at  -2900  km  from  the  source  provides  a  sufficient  coda 
time  window  before  the  onset  of  the  S-wave  phase. 

Crustal  and  coda  attenuation.  It  is  important  to  differentiate  between  two  different  Q  values  used  in  this  study. 

The  intrinsic  or  elastic  crustal  attenuation  describes  the  loss  of  seismic  energy  in  the  crust  to  internal  friction,  such  as 
damping  by  fluids  in  cracks  and  pores  or  friction  along  grain  boundaries.  It  is  denoted  by  a  Q  value,  where 

„  2  kE 

Q  — - ,  with  E  being  the  energy  and  AE  -  the  energy  dissipation  during  one  wave  cycle  (Aki  and  Richards, 
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2002).  By  contrast,  0coda  measures  the  coda  amplitude  decay  rate,  determined  from  the  decay  slope  of  amplitude  A 
versus  time  t  from  a  semi-log  plot  of  the  seismic  trace  envelope  in  Eq.  (2).  By  taking  the  natural  logarithm: 

i  j  i  j  7tf  -7tf 

111  A  —  in  Aq - 1 ,  we  obtain  a  straight  line  with  slope  S  — -  in  the  ( time,\nA )  domain.  Thus,  the 

Qcoda  Qcoda 


frequency-dependent  Qcoda(f)  is  defined  as: 


,  where  slope  s  is  measured  from  the  logarithm  of  the 


seismic  trace  envelope. 


Numerical  Coda  Model.  The  scattering  model  represents  crustal  heterogeneities  as  scatterers  distributed  within 
Earth’s  crust.  The  resultant  scattered  coda  intensity  U recorded  at  a  receiver  at  time  t  is  an  integral  over  volume  V 
containing  all  scatterers  (Figure  1): 


U{r,t)  =  \dts\\d2rs'¥(rs)U: 

s 


(rs,ts)G(r-rs,t-ts). 


(4) 


Figure  1.  Surface  scattering  model  (Eq.  4).  Seismic  energy  originates  at  the  source,  scatters  from  randomly 
and  uniformly  distributed  surface  points,  and  is  detected  at  receiver.  The  source  and  Green’s 
functions  are  simulated  by  1-D  reflectivity  modeling. 


For  a  given  coda  time  t  and  velocity  of  a  particular  scattered  mode  (e.g.,  v  =  2.9  -  3.5  km/s  for  a  scattered  Lg ), 
contributions  of  scattered  energy  originate  from  an  elliptical  ring  surrounding  the  receiver  (Morozov  and  Smithson 
2000).  The  scattering  area  within  the  ring  increases  with  time  due  to  its  increasing  radius  and  partly  compensates  the 
energy  decay  from  the  geometrical  spreading.  For  example,  for  a  plane  of  constant  scattering  potential  and  no 
intrinsic  attenuation,  coda  energy  would  not  decay,  provided  only  surface-wave  modes  contribute  to  the  scattered 
Green’s  function. 

Vp  [km/s]  Vp  [km/s] 

456789  456789 


Figure  2.  Upper  250  km  of  1-D  VP  velocity  models  used  in  creation  of  synthetic  seismograms,  a)  Models  based 
on  PNE  Quartz  results  by  Morozova  et  al.  (1999):  Quartz  (solid  black  line),  Warm  Quartz  (dotted 
black),  and  Complex  Crust  (gray);  b)  iasp91  model. 

Synthetic  seismic  sections  created  using  the  reflectivity  method  (Fuchs  and  Muller,  1971)  were  used  as  source  and 
Green’s  functions  in  several  1-D  crustal  and  upper  mantle  models  (Figure  2).  The  original  FORTRAN  program  by 
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K.  J.  Sandmeier  was  modified  to  handle  larger  computations,  parallelized,  and  incorporated  into  our  data  processing 
system  (SI A;  Chubak  and  Morozov,  2006)  allowing  seamless  filtering,  inverting,  saving,  and  plotting  the  results. 

A  detailed  description  of  this  tool  can  be  obtained  from  the  SI  A  web  service: 
http://seisweb.usask.ca/SIA/ps.php?doc=reflect. 

Computations  were  performed  on  a  66-processor  1.6-GHz  Opteron  cluster  and  produced  800-second,  3-component 
synthetic  records  sampled  at  20-ms  intervals  and  output  at  10-km  intervals  from  near-0  to  3500  km  from  the  source. 
Modeling  frequency  band  was  0.2-20  Hz  using  a  “delta- function”  source  function  suitable  for  spectral 
measurements.  A  sufficiently  dense  phase  velocity  spectrum  was  selected  to  avoid  frequency  aliasing  during 
numerical  mode  summations.  As  with  any  implementations  of  the  propagator  matrix  method  (Aki  and  Richards, 
2002),  all  P/SV  mode  conversions  and  multiples  were  accounted  for  in  the  modeling. 

We  used  four  different  velocity  models  to  compute  the  synthetics  (Figure  2),  including  the  global  iasp91  model 
(Kennett  and  Engdahl,  1991),  consisting  of  a  simple  three-layer  crust  and  a  mantle  without  strong  gradients  and  low- 
velocity  zones.  However,  the  iasp91  model  is  too  slow  for  platform  areas  of  northern  Eurasia,  and  the  remaining 
three  models  were  based  on  the  detailed  models  derived  from  the  NE  profile  Quartz  and  considered  the  best  match 
for  the  study  area  (Mechie  et  al.,  1993;  Morozova  et.  al.,  1999).  The  “Quartz”  model  contains  a  3-gradient  layer 
crust  overlain  by  a  3-km-thick  sediment  layer  and  complex  mantle  with  low  velocity  zones  at  1 10-  and  2 10-km 
depths  (Figure  2).  Another  “Complex  Crust”  model  has  a  high-contrast,  5 -layer  crust  with  a  somewhat  exaggerated 
attenuative  layer  ( Qs  =  10)  on  its  top  and  the  Quartz  mantle  (Figure  2b).  Model  “Warm  Quartz”  (Figure  2a)  was  an 
attempt  to  simulate  a  high-heat  flow  regime  in  the  Quartz  model  by  applying  a  temperature-related 
negative  velocity  gradient  within  the  crust  (Christensen  and  Mooney,  1995).  Note  that  all  three  Quartz- 
based  models  show  more  complex  waveforms  (including  crustal  Lg  phases  in  addition  to  the  P  and  S 
phases)  compared  with  the  iasp91  (Figure  2).  In  order  to  collect  data  on  the  sensitivity  of  the  results  to 
crustal  velocities  alone,  we  also  repeated  the  Quartz  simulations  after  applying  ±5%  perturbations  of 
crustal  velocities. 

Multiple  simulations  were  performed  for  each  model,  with  P-  and  5-wave  quality  factors  for  layers  in  the 
crustal  part  of  the  model  set  equal — QP  =  2 Q  and  Qs  =  Q ,  respectively,  with  values  of  Q  ranging  from  100 
to  1000  to  encompass  the  expected  Earth’s  whole-crust  Qs  values.  This  range  of  Qs  corresponds  to  Lg  Q 
values  of  -300-1000  measured  from  19  Russian  PNEs  by  Li  et  al.  (submitted  to  the  Bulletin  of  the 
Seismological  Society  of  America  [BSSA])  and  extended  to  low-0  values  to  investigate  crustal  effects. 
Mantle  attenuation  in  each  velocity  model  was  fixed  using  the  mantle  Qp  values  derived  from  P-wave 
Quartz  PNE  data  (Morozov  et  al.,  1998b).  The  teleseismic  P  coda  was  measured  at  2900-km  offsets  to 
allow  a  significant  time  interval  in  which  to  perform  the  measurement  before  the  onset  of  the  Pg  and  S 
waves.  Numerical  evaluation  of  integral  (4)  was  performed  by  using  Monte-Carlo  sampling  in  order  to 
avoid  spurious  coherency  that  could  be  caused  by  spatial  aliasing  when  using  a  regular  integration  grid. 

As  above,  the  synthetic  coda  generation  was  implemented  in  an  SIA  tool  sharing  much  of  its  libraries 
with  other  codes  (http  ://seis web .usask. ca/SI A/ps .php?doc=syncoda) . 

Using  the  resulting  coda  records,  log-amplitude  slope  measurements  were  performed  on  coda  envelopes 
in  the  manner  of  Morozov  and  Smithson  (2000)  (Figure  3).  Both  L2-  (root  mean  square  [RMS])  and  Lr 
norm  fitting  (e.g.,  Aster  et  al.,  2002)  were  tried,  with  uncertainties  measured  using  50%  jack-knife 
resampling.  Owing  to  large  sampling  volumes  and  stable  coda  envelopes,  both  the  differences  between 
the  L2  and  Li  norms  and  the  corresponding  slope  uncertainties  were  negligible  for  the  subsequent 
analysis. 

Modeling  Results.  We  performed  measurements  of  the  teleseismic  P-wave  codas  for  each  of  the  four 
velocity  models.  Coda  traces  were  band-pass  filtered  to  center  frequencies  of  0.4,  0.79,  1.58,  3.16,  4.48, 
6.32,  8.94,  12.6,  and  17.8  Hz  to  examine  the  log-amplitude  slopes  as  functions  of  frequency.  Figure  4 
shows  the  results  for  models  Quartz  and  Complex  Crust  with  Qs  =  Q  =  500. 

In  Figure  4,  coda  slope  values  show  clear  linear  dependencies  to  about  7  to  10  Hz,  after  which  the  fidelity 
of  modeling  appears  to  drop,  and  the  slopes  flatten  out  (attenuation  is  replaced  by  frequency-independent 
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numerical  noise  of  unclear  nature).  By  fitting  straight  lines  to  the  points  within  the  0  to  7  Hz  range,  we 
find  that  the  intercepts  at /=  0  Hz  are  non-zero,  y~  0.8*  10"2,  and  y~  1.6-10"2  for  the  two  models  and 
measure  Qco da  from  their  slopes  (see  Eq.  3)  as  QCOda  ~  800  and  600,  respectively,  Note  that  yis  clearly 
positive  in  all  cases  considered  here,  indicating  that  the  geometrical  spreading  is  only  partly  compensated 
by  surface  scattering.  Also  note  the  flattening  of  slopeif)  dependencies  at  higher  frequencies  due  to  the 
reduction  in  the  signal/noise  ratio  and  outliers  at  lower  frequencies,  most  likely  associated  with  tuning 
effects  within  the  crustal  layers  (Figure  4). 


200  400 


600  800 
Time  [s] 


Figure  3.  Log-amplitude  synthetic  coda  for  Quartz  model  at  selected  frequencies  (labeled).  Straight 
lines  indicate  the  measured  amplitude  fits.  Note  the  progressive  change  in  coda  shapes  with 
increasing  frequencies. 
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Figure  4.  Measured  log-amplitude  coda  slopes  from  models  Quartz  and  Complex  Crust  with  crustal 
Qs  =  500.  Note  the  linearly  decreasing  (steepening)  coda  slopes  with  frequencies  increasing 
to  ~7-10  Hz.  Solid  gray  lines  indicate  linear  regressions  using  formula  (3):  y=  0.8*10"2  and 
£?coda  =  800  for  Quartz,  and  y=  1.6-10"2  and  Qc oda  =  600  for  Complex  Crust  models. 
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Alternatively,  when  converted  to  Qco da  by  using  Eq.  (2),  the  values  of  slope  result  in  frequency-dependent 
2coda(/)  (Figure  5),  which  one  could  also  fit  using  a  dependence  of  type  (1)  with  Q0  ~  250  at  1  Hz  and  77  ~ 
0.5  [see  the  discussion  of  conversion  to  (77,  Qo )  below].  Although  this  interpretation  could  be  acceptable 
within  the  synthetic  or  real  data  uncertainties,  our  model  contains  no  mechanisms  to  support  an  in  situ 
frequency-dependent  attenuation.  Therefore  the  observed  QC0da(f)  (Figure  5)  should  be  viewed  as  an 
“apparent”  quantity  or  as  an  observational  artifact. 

By  selecting  linear  In (A(t))  dependencies  (3)  for  each  frequency,  we  found  that  the  values  of  ^stayed 
approximately  constant  in  all  cases  (values  shown  in  labels  in  Figure  6).  These  observations  were 
expected  (see  Discussion ),  as  y  should  principally  depend  on  the  crustal  structure.  By  using  these 
constant  afs  and  for  any  frequency  (Figure  6),  the  observed  coda  slopes  can  be  transformed  into  the 
following  ratios: 


a{<2)  = 


Qcoda 

Q 


4 


-Q 


d\nA(t,f ) 


dt 


+  y 


(5) 


Note  that  the  resulting  QcodJQ  ratio  is  nearly  frequency-independent,  greater  than  1,  and  approaches  1  for 
high  Q.  Because  of  these  properties,  o(Q)  can  be  tabulated  and  approximated  from  numerical  modeling 
(lines  in  Figure  6)  and  utilized  as  a  “calibration  curve”  for  inverting  the  Qco da  observations  for  crustal  Q. 
Interestingly,  functions  are  very  close  (within  -0.2  difference)  for  all  three  Quartz-based  models,  despite 
the  strong  difference  in  y  for  the  Complex  Crust  model  (Figure  6).  However,  the  a(Q)  for  iasp91  model 
is  strongly  different  from  those  for  Quartz  in  its  shape,  values,  and  stability  in  respect  to  frequency 
variations  (Figure  6). 


Frequency  [Hz] 


Figure  5.  Measured  log-amplitude  coda  slopes  transformed  to  apparent  coda  Qcoaa(f)  according  to 
relation  (2).  Note  the  strong  frequency  dependence  (//M).5)  of  QCOda(f)  despite  the  constant- 
Q  modeled  crustal  rheology. 
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Finally,  both  ^and  the  functional  form  of  a(Q)  should  depend  on  the  velocity  and  density  structures  of 
the  models.  We  tested  this  dependence  by  using  four  different  types  of  models  (Figure  2)  and  also  by 
adding  velocity  perturbations  to  Quartz  model.  The  results  suggest  that  ^and  a(Q )  indeed  correlate  with 
the  structural  model  styles  and  crustal  velocities,  as  will  be  further  elaborated  on  in  Discussion. 
Mapping  the  (j^Qcoda)  parameterization  to  ( Q{ 0,  rj).  In  most  practical  cases,  the  time-frequency 
amplitude  decays  A(f,  t)  can  be  satisfactorily  matched  in  either  (y  Qco da)  or  (Q0,  if)  forms,  and 
approximate  mapping  between  them  can  be  developed.  Consider  an  attenuation  process  (e.g.,  coda) 
governed  by  Eq.  (3)  but  parameterized  by  Eqs.  (2)  and  (1).  Therefore,  for  the  time-dependent  part  of 
lnA(f, ’ t ),  we  have 


where  q  is  a  constant.  For  a  specific  range  of  frequencies  and  coda  times,  transformation  (y  Qcoda)  — >  ( Qo , 
rj)  can  be  found  numerically  by  least-squares  fitting  the  logarithm  of  Eq.  (6): 


lng  +  (l-77)ln /  ~  In 


y+ 

v  a 


¥ 


coda  J 


(7) 


and  solving  it  for  \nq  and  (I-77).  The  resulting  (Q0,  if)  parameters  for  our  numerical  experiments  are 
shown  in  Figure  7.  Note  that  77  is  principally  controlled  by  the  yQcoda product,  with  Q0  =  Qco da  when  y=  0 
and  decreasing  with  increasing  y  (Figure  7).  Thus,  high  77  ~1  should  be  expected  in  areas  of  low 
attenuation,  with  high  geometric  spreading,  or  when  measurements  are  carried  out  at  lower  frequencies. 
Therefore,  stable  platforms  (such  as  the  Siberian  Craton  in  eastern  Russia)  could  exhibit  high  apparent 
frequency  dependence  of  Qc 0da  simply  on  account  of  their  high  Q.  More  detail  on  this  can  be  found  in 
Morozov  et  al.  (submitted  to  BSSA). 


QS=Q  Q,  =  Q 

Figure  6.  a{Qs)  from  Eq.  (5)  for  the  four  models  of  this  study  and  different  frequencies  (labeled 
symbols),  obtained  directly  from  coda  slopes  by  using  ^values  in  Eq.  (3)  shown  in  the 
labels.  The  lines  represent  the  suggested  empirical  a{Qs)  dependences  for  each  velocity 
structure,  picked  by  hand. 
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PNE  data  example.  The  use  of  (y  Qco da)  description  allows  us  to  resolve  the  following,  recent,  and 
hitherto  puzzling  observation  from  the  PNE  profiles.  Within  the  east  European  platform,  measurements 
from  Quartz  profile  resulted  in  Qc oda(~2  Hz)  =  380  and  0coda(~5  Hz)  ~  430,  which  was  interpreted  in  the 
power-law  form  (1)  as  Qcoda(f)  ~  210  f 03  (Morozov  and  Smithson,  2000).  However,  within  the  Siberian 
Craton  (PNE  Kimberlite-3),  practically  frequency-independent  coda  amplitude  decays  were  observed, 
corresponding  to  QCOda(f)  ~  1050/1 0  (Morozov  et  al.,  2006).  Although  the  attenuation  was  expected  to  be 
low  within  the  Siberian  Craton,  its  strong  frequency  dependence  is  vastly  different  from  the  east  European 
platform  and  is  surprisingly  high  even  for  tectonically  active  areas  (Erickson  et  al.,  2004). 


The  explanation  of  the  above  in  the  (  %  Qco da)  form  is  quite  simple.  From  the  Quartz-4  data,  the  logarithms 
of  coda  amplitudes  can  be  interpreted  by  using  y  ~  0.75T0"2  s"1  and  frequency-independent  Qco da  *  850, 
with  an  estimated  range  of  uncertainty  of  780  to  960  (Figure  8).  From  Kimberlite-3  records,  j^turns  out  to 
be  nearly  the  same,  and  gcoda  ~  2500±300  (Figure  8).  Thus,  the  Siberian  Craton  should  occupy  a  region  in 
the  upper  parts  of  the  diagrams  in  Figure  7,  where  rj~  1.  Note  that  both  values  of  j'also  agree  remarkably 
well  with  the  numerical  simulations  based  on  Quartz  velocity  model  data  above  (Figure  4).  The  values  of 
Q  within  the  Siberian  Craton  are  very  high,  in  agreement  with  observations  of  short-period  Pg  waves 
propagating  to  over  1600  km  from  the  PNEs  (Morozov  et  al.,  2006).  These  are  likely  the  longest- 
propagating  short-period  Pg  observed  anywhere  on  Earth. 
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Figure  7.  Mapping  of  (y  Qco da)  to  ( Q0 ,  rj)  approximated  for /requencies  of  1  to  10  Hz.  The 

hyperbolas  drawn  in  both  plots  correspond  to  QCOd*Y=  tpIHz.  Note  that  ^increases  with 
increasing  Qcoday 

CONCLUSIONS  AND  RECOMMENDATIONS 

Frequency-dependent  attenuation  often  observed  in  regional  coda  measurements  can  be  interpreted  in  two 
ways:  by  frequency-dependent  coda  QC0da(f)  or  by  means  of  effective,  uncompensated  geometrical 
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spreading  (y)  and  frequency-independent  coda  attenuation  (Qco  da).  The  second  of  these  alternatives  offers 
a  number  of  practical  advantages:  1)  it  is  better  justified  theoretically  and  does  not  require  unusual  crustal 
rheologies  or  inverted  scaling  properties;  2)  its  parameters  can  be  quantitatively  related  to  the  in  situ 
velocity/density  (y)  and  attenuation  (Qco da)  crustal  structures;  3)  the  resulting  Qcoda  appears  to  correspond 
to  crustal  Qs  and  can  be  directly  correlated  to  the  observed  Lg  Q\  4)  the  frequency-independent 
parameters  allow  easier  comparisons  between  different  areas  and  studies;  and  5)  as  stable  characteristics, 
these  parameters  are  likely  to  improve  coda  magnitude  analysis  and  coda  Q  regionalization. 

Frequency-dependent  inelasticity  due  to  spectral  variations  of  scatterer  density  is  certainly  possible  and 
likely  to  be  present  in  Earth’s  crust.  However,  the  objective  of  this  paper  was  to  show  that  it  may  not  be 
as  pervasive  as  suggested  by  many  frequency-dependent  Q  observations.  Within  the  available  data 
quality,  observations  of  apparent  Q(f)  at  short  periods  could  often  be  successfully  explained  by 
combinations  of  geometrical  spreading  (or  turbidity  in  elastic  scattering — Dainty,  1981)  and  frequency- 
independent  attenuation.  Separation  of  these  physical  parameters  of  the  crust  provides  several  insights  for 
interpretation  and  allows  linking  the  seismic  velocity,  structure,  coda,  and  Lg  Q  results.  This  conclusion  is 
almost  certainly  the  case  with  PNE  data,  and  it  may  also  apply  to  a  variety  of  other  seismic  datasets  (e.g., 
surface  waves). 

Numerical  modeling  of  PNE  coda  wavefields  argues  strongly  in  favor  of  the  new  (y  Qco da)  interpretation 
and  illustrates  the  mechanisms  leading  to  such  an  apparent  Q(f).  The  inverted  Qco da  values  can  be 
systematically  mapped  into  the  true  S- wave  attenuation  factors  within  the  crust.  By  contrast,  when 
interpreted  using  the  traditional  approach,  the  synthetic  coda  shows  strong  spurious  frequency-dependent 

2coda* 

PNEs  in  Russia  in  two  different  areas  were  found  to  have  similar  ^values,  although  strongly  different 
2coda-  We  also  suggest  that  parameter  j'could  provide  a  better,  stable,  and  transportable  discriminant  for 
differentiating  between  tectonic  types  of  the  crust.  From  a  limited  preliminary  example,  this  discriminant 
may  relate  to  the  presence  of  localized  attenuative  or  scattering  structures  within  the  crust  and  be 
applicable  to  both  North  America  and  Northern  Eurasia. 


Figure  8.  Log-amplitude  coda  slopes  in  real  data  from  PNEs  Quartz-4  (recorded  within  the  east 
European  platform,  at  ~2600-km  source-receiver  offset;  same  data  as  in  Morozov  et  al., 
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2002)  and  Kimberlite-3  (recording  within  the  Siberian  Craton  at  ~1100-km  offsets).  Note 
the  different  Qc oda  values  (lines  and  labels),  whereas  the  ^values  are  similar  and  close  to 
those  modelled.  Also  note  that  because  of  lower-magnitude  source  and  longer  recording 
distance,  the  usable  frequency  band  is  narrower  for  Quartz  PNE  records. 
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